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Abstract 

QCD sum rules for the magnetic moments of A ++ and £l~ are derived using 
the external field method. They are analyzed by a Monte-Carlo based pro- 
cedure, using realistic estimates of the QCD input parameters. The results 
are consistent with the measured values, despite relatively large errors that 
can be attributed mostly to the poorly-known vacuum susceptibility %■ K is 
shown that a 30% level accuracy can be achieved in the derived sum rules, 
provided the QCD input parameters are improved to the 10% level. 
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The QCD sum rule method is a powerful tool in revealing the deep connection between 
hadron phenomenology and QCD vacuum structure via a few condensate parameters. The 
method has been successfully applied to a variety of problems to gain a field-theoretical 
understanding into the structure of hadrons. Calculations of the nucleon magnetic moments 
in the approach were first carried out in Refs. [0] and [^]. They were later refined and 
extended to the entire baryon octet in Refs. . On the other hand, the magnetic moments 
of decuplet baryons were less well studied within the same approach. There were previous, 
unpublished reports in Ref. [§] on A ++ and Q~ magnetic moments. The magnetic form factor 
of A ++ in the low Q 2 region was calculated based on a rather different technique ||. In 



recent years, the magnetic moment of £7 has been measured with remarkable accuracy []T0 



fj,Q- = —2.02 ± 0.05 /i at. The magnetic moment of A ++ has also been extracted from pion 



bremsstrahlung ||11|| : /xa++ — 4.5 ± 1.0 /ijv- The experimental information provides new 
incentives for theoretical scrutiny of these observables. 

In this letter, we present an independent calculation of the magnetic moments of A ++ 
and Q~ in the QCD sum rule approach. The goal is two-fold. First, we want to find out 
the applicability of the method as an alternative way of understanding the measured /ia++ 
and fiQ- , and hope to gain some insights into the internal structure of these baryons from a 
nonperturbative-QCD perspective. Second, we want to achieve some realistic understanding 
of the uncertainties involved in such a determination by employing a Monte-Carlo based 
analysis procedure. This will help find possible ways for improvement. The calculation is 
algebraically more involved than the octet case since one has to deal with the more complex 
spin structure of spin-3/2 particles, but presents no conceptual difficulties. 

Consider the two-point correlation function in the QCD vacuum in the presence of a 
constant background electromagnetic field F^ u : 

n^(p) = i J d'x e tp - x (0 | T{ Va (x) %(0) } | 0) F , (1) 
where i] a is the interpolating field for the propagating baryon. The subscript F means 
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that the correlation function is to be evaluated with an electromagnetic interaction term: 
Cj = —A^J 11 , added to the QCD Lagrangian. Here is the external electromagnetic 
potential and is given by A^ix) = — \F llv x u in the fixed-point gauge, and = e q (pi^q the 
quark electromagnetic current. The magnetic moments can be obtained by considering the 
linear response of the correlation function to the external field. The action of the external 
field is two-fold: it couples directly to the quarks in the baryon interpolating fields, and it 
also polarizes the QCD vacuum. The latter can be described by the introduction of vacuum 
susceptibilities. 

The interpolating field is constructed from quark fields, and has the quantum numbers 
of the baryon in question. For A ++ and Q~ , they are given by 

n* ++ {x) = e abc (u aT {x)C la u\x)) u%x), (2) 

rfc-(x) = e abc (s aT (x)C la s b (x)) s c (x), (3) 

where C is the charge conjugation operator, and the superscript T means transpose. The 
ability of a interpolating field to annihilate the ground state baryon into the QCD vacuum 
is described by a phenomenological parameter X B (called current coupling or pole residue), 
defined by the overlap 

(0 \n a | Bps) = XsUaip, s), (4) 

where u a is the Rarita-Schwinger spin-vector. 

The QCD sum rules are derived by calculating the correlator in ([[]) using Operator- 
Product-Expansion (OPE), on the one hand, and matching it to a phenomenological repre- 
sentation, on the other. The calculation is similar in spirit to the octet case (see Ref. [£§, 
for example), but with the added complexity of spin-3/2 structures. A direct evaluation 
led to numerous tensor structures, not all of them independent. The dependencies can be 
removed by ordering the gamma matrices in a specific order (we chose p^alnlvlp where the 
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hat notation means p = £> M 7 M ). After a lengthy calculation, a sum rule involving only the 
magnetic moment can be isolated at one of the structures: ^aF^a^pp. Here we give the 
final results, for A ++ : 

-e^E.M 4 + -e u bL m ~ -e uX a 2 L 12 ' 27 - -e u K v a 2 L 28 / 27 -L 
28 56 7 7 M 2 

= ' x i(w + A ) e ' M ' B,M '- (5) 

and for 

A esL 4/ 27 £ iM 4 _ ^ esf(f)msXaL -^ EoM 2 + l es bL 4/27 ~ fe s fm s aL^ 
~e s M2K + Z)m.aL"™ - \e s f 2 <\> X a 2 L^ - ^e s f 2 n v a 2 L 2 ^ ^ 

-y^/ 2 ^ + ^ L28/27 w 2 + \^ 2 ^< a2L ~ 2l21 w 

_ A es/msm 2 aL -10/27 1 + 1 /2 2 2 L 14/27 J_ 

28 M 2 12 M 4 

= A|(^ + A)e--I/- 2 . (6) 

In these equations, the magnetic moment /i# is given in particle's natural magnetons. The 
various symbols are defined as follows. The condensates are represented by a = —(2n) 2 (uu), 
b = (g 2 G 2 ), (ug c a ■ Gu) = —ml (uu), and the coupling \ B = (27t) 2 Ab. The factors e u = 2/3 
and e s = —1/3 are quark charges in units of electric charge. The vacuum susceptibil- 
ities x, k and f are defined by (qcr^q) F = e q x(qq)F^ v , (qg c G^q) F = e q n(qq)F^ v , and 
(qg^^p\G px ^q) F = ieqti^F^. The parameters / = (ss)/ (uu) = (sg c a ■ Gs)/ (ug c a ■ Gu) 
and = Xs/x — k s /k = £ s /£ account for flavor symmetry breaking of the strange quark. 
Possible violation of the four-quark condensate is considered by the parameter k v as defined 
in (uuuu) = k v (uu) 2 . The anomalous dimension corrections of the various operators are 
taken into account via the factor L = \a s (^ 2 ) / a s (M 2 )} = hi(M 2 / k 2 QCD ) / h\(n 2 / h 2 QCD ) , 
where \x = 500 MeV is the renormalization scale and Aq C d is the QCD scale parameter. 
As usual, the pure excited state contributions are modeled using terms on the OPE side 
surviving M 2 — > oo under the assumption of duality, and are represented by the factors 



E n (x) = 1 — e~~ x J2 n x n /n\ with x = w 2 jM\ and w an effective continuum threshold. The 
parameter A accounts for all contributions from the transitions caused by the external field 
between the ground state and the excited states. Such contributions are not exponentially 
suppressed relative to the ground state double pole and must be included. The presence of 
such contributions is a general feature of the external field technique. 

Let us note in passing that since Q~ and A ++ are simply related by the interchange of 
quark flavors u <-> s, one can verify that sum rule (|^) reduces to sum rule ((^) if one sets 
e s — > e u , m s = 0, f = 1, and = 1. 

To analyze the sum rules, we use a Monte-Carlo based procedure first developed in 
Ref. [ pfl . The basic steps are as follows. First, the uncertainties in the QCD input parame- 
ters are assigned. Then, randomly-selected, Gaussianly-distributed sets are generated, from 
which an uncertainty distribution in the OPE can be constructed. Next, a \ 2 minimization 
is applied to the sum rule by adjusting the phenomenological fit parameters. This is done 
for each QCD parameter set, resulting in distributions for phenomenological fit parameters, 
from which errors are derived. Usually, 100 such configurations are sufficient for getting 
stable results. We generally select 1000 sets which help resolve more subtle correlations 
among the QCD parameters and the phenomenological fit parameters. 

The Borel window over which the two sides of a sum rule are matched is determined 
by the following two criteria. First, OPE convergence: the highest-dimension-operators 
contribute no more than 10% to the QCD side. Second, ground-state dominance: excited 
state contributions should not exceed more than 50% of the phenomenological side. The 
first criterion effectively establishes a lower limit, the second an upper limit. Those sum 
rules which do not have a valid Borel window under these criteria are considered unreliable 
and therefore discarded. 



The QCD input parameters and their uncertainty assignments are given as follows |L3 . 



The condensates are taken as a = 0.52±0.05 GeV 3 , b = 1.2±0.6 GeV 4 , and = 0.72±0.08 
GeV 2 . For the factorization violation parameter, we use k v = 2 ± 1 and 1 < k v < 4. The 



QCD scale parameter is restricted to Aq(7£>=0.i5±0.04 GeV. The vacuum susceptibilities 
have been estimated in studies of nucleon magnetic moments [@-[§], but the values vary in 
a wide range depending on the method used. Here we take their median values with 50% 
uncertainties: x = -6.0 ± 3.0 GeV~ 2 , « = 0.75 ± 0.38, and f = -1.5 ± 0.75. Note that \ is 
almost an order of magnitude larger than k and £, and is the most important of the three. 
The strange quark parameters are placed at m s = 0.15 ± 0.02 GeV, / = 0.83 ± 0.05 and 
= 0.60 ± 0.05 [f|. These uncertainties are assigned conservatively and in accord with the 
state-of-the-art in the literature. While some may argue that some values are better known, 
others may find that the errors are underestimated. In any event, one will learn how the 
uncertainties in the QCD parameters are mapped into uncertainties in the phenomenological 
fit parameters. 

To illustrate how well a sum rule works, we first cast it into the subtracted form, lis = 
A^/i_Be~ M s/ M2 ; then plot the logarithm of the absolute value of the two sides against the 
inverse of M 2 . In this way, the right-hand side will appear as a straight line whose slope is 
— M B and whose intercept with the y-axis gives some measure of the coupling strength and 
the magnetic moment. The linearity (or deviation from it) of the left-hand side gives an 
indication of OPE convergence, and information on the continuum model and the transitions. 

To extract the magnetic moments, a two-stage fit was performed. First, the correspond- 
ing chiral-odd mass sum rule, as obtained previously in Ref. |]I3"|| , was fitted to get the mass 
Mb, the coupling \ B and the continuum threshold w±. Then, Mb and \ 2 B were used in the 
magnetic moment sum rule for a three-parameter fit: the transition strength A, the contin- 
uum threshold and the magnetic moment /xg. Note that w± and W2 are not necessarily 
the same. We impose a physical constraint on both w± and W2 requiring that they are larger 
than the mass, and discard QCD parameter sets that do not satisfy this condition. The 
above procedure is repeated for each QCD parameter set until a certain number of sets are 
reached. In the actual analysis of sum rules (|5|) and (|6]), however, we found that such a 
full search was unsuccessful: the search algorithm consistently returned W2 either zero or 
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smaller than Mb- This signals insufficient information in the OPE to completely resolve the 
spectral parameters. To proceed, we fixed W2 at w±, which is a reasonable and commonly 
adopted choice in the literature, and searched for A and /xg. 

Fig. |l| shows the match for the sum rules fl5|) and (H). The extracted results are given 
in Table |. Relatively large errors in the magnetic moments are found, approaching 100%. 

~~ 1 1 1 1 1 1 1 1 

T T 

A ++ 



j 1 1 1 1 

0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

M~ z (GeV~ 2 ) 

FIG. 1. Monte-Carlo fits of the magnetic moment sum rules @ and (^). Each sum rule is 
searched independently. The solid line corresponds to the ground state contribution, the dotted 
line the rest of the contributions (OPE minus continuum minus transition). The error bars are 
only shown at the two ends for clarity. 

But the sign and order of magnitudes are unambiguous when compared to the measured 
values. The situation is consistent with a previous finding on qa 0] regarding three-point 
functions. It is interesting to observe that the distribution of errors is not uniform throughout 
the Borel window, with the largest errors at the lower end where the power corrections are 
expected to become more important. The quality of the match deteriorates for Q~ in this 
region, signaling probably insufficient convergence of the OPE. The match for A ++ is good 
in the entire Borel region, despite the large uncertainties. To gain some idea on how the 
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TABLE I. Monte-Carlo analysis of the QCD sum rules for the magnetic moments of A ++ and 
il~ . The third column represents the percentage contribution of the excited states and transitions 
to the phenomenological side at the lower end of the Borel region (it increases to 50% at the upper 
end). The second row for each sum rule corresponds to results with reduced, uniform 10% errors 
assigned to all the QCD input parameters. The uncertainties were obtained from consideration of 
1000 QCD parameter sets. In the event that the resultant distribution is not Gaussian, the median 
and asymmetric deviations are reported. 



Sum Rule 


Region 
(GeV) 


Cont 

(%) 


w 
(GeV) 


A 

(GeV" 2 ) 


MB 

(Pn) 


d) for A++ 

d) for n- 


0.765 to 1.47 
0.765 to 1.47 
0.747 to 1.66 
0.747 to 1.66 


9 
8 
7 
6 


1.65 
1.65 
2.30 
2.30 


0.53 ±0.77 
0.35 ±0.37 
-0.15 ±0.14 
-0.10 ±0.05 


3.60±H 

4 1O+1.40 

-1 25+ 112 
1 4 q+0.40 

-i.4y_ 049 


uncertainties chan 


.ge with the input, 


we also analyzed the 


sum rules by adjusting 


; the error 



estimates individually. We found large sensitivities to the susceptibility x- 111 f&ct, most 
of the errors came from the uncertainty in \. We also tried with reduced error estimates 
on all the QCD input parameters: 10% relative errors uniformly. The results are given 
in Table | as a second entry. It leads to about 30% accuracy on the magnetic moments. 
Further improvement of the accuracy by reducing the input errors is beyond the capability 
of these sum rules as it will lead to unacceptably large x 2 I^df [0- For that purpose, one 
would have to resort to finding sum rules that depend less critically on \ and have better 
convergence properties. 

A comparison with those from other calculations and the experimental data is compiled 
in Table |||. The results with 10% errors from the QCD sum rule method are used in the 
comparison. They are consistent with data within errors, although the central value for Q~ 
is somewhat underestimated. The result for A ++ is consistent with other calculations, while 
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TABLE II. Comparisons of magnetic moments from various calculations: this work (QCDSR), 
lattice QCD (Latt) [fig] , chiral perturbation theory (%PT) |]l6f , light-cone relativistic con- 
stituent quark model (RQM) pTfl , simple non-relativistic constituent quark model (NQM), chiral 
quark-soliton model (xQSM) ]i~8|] . All results are expressed in units of nuclear magnetons. 



Baryon 


Exp. 


QCDSR 


Latt 




RQM 


NQM 


XQSM 


A++ 


4.5 ± 1.0 


4.13 ± 1.30 


4.91 ± 0.61 


4.0 ± 0.4 


4.76 


5.56 


4.73 


n~ 


-2.024 ± 0.056 


-1.49 ± 0.45 


-1.40 ± 0.10 




-2.48 


-1.84 


-2.27 



for ft it is closer to lattice QCD calculations. 

In conclusion, we have demonstrated that the magnetic moments of A ++ and Q~ can be 
understood from the QCD sum rule approach, despite large errors that can be traced to the 
uncertainties the QCD input parameters. A 30% accuracy can be achieved in the derived 
sum rules with improved estimates of the QCD input parameters, preferably on the 10% 
accuracy level. In particular, large sensitivities to the vacuum susceptibility x are found. 
Better estimate of this parameter is needed. Extension of the calculations to other members 



of the decuplet appears straightforward, and is under way |19| . There we hope to address 
the issues involved in greater detail. 

It is a pleasure to thank D.B. Leinweber for providing an original version of his Monte- 
Carlo analysis program and for helpful discussions. This work was supported in part by 
U.S. DOE under Grant DE-FG03-93DR-40774. 
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